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Abstract. 

In this paper the sensitivity of a future kilometer-sized neutrino detector to detect 
and measure the diffuse flux of high energy neutrinos is evaluated. Event rates 
in established detection channels, such as muon events from charged current 
interactions or cascade events from v & and v r interaction, are calculated using a detailed 
Monte Carlo simulation. Neutrino fluxes as expected from prompt charm decay in the 
atmosphere or from astrophysical sources such as Active Galactic Nuclei (AGN) are 
modeled assuming power laws. The ability to measure the normalization and slope of 
these spectra is then analyzed. 

It is found that the cascade channel generally has a high sensitivity for the detection 
and characterization of the diffuse flux, when compared to what is expected for the 
upgoing- and downgoing-muon channels. A flux at the level of the Waxman-Bahcall 
upper bound should be detectable in all channels separately while a combination of 
the information of the different channels will allow detection of a flux more than one 
order of magnitude lower. Neutrinos from the prompt decay of charmed mesons in 
the atmosphere should be detectable in future measurements for all but the lowest 
predictions. 



PACS numbers: 14.60.Pqx,96.40.Tv,95.85.Ry 
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1. Introduction 

High energy extra-terrestrial neutrinos have so far escaped their detection, and there is 
a considerable effort invested to change this situation soon pQ. IceCube, a Cherenkov 
detector with a volume of a cubic kilometer is currently being installed at the South 
Pole [2] and plans for a northern hemisphere detector of similar size are maturing 
Such detectors are expected to measure high energy neutrinos [E v > 100 GeV) from 
various sources of which only a few might be resolved by directional information. The 
unresolved diffuse (or isotropic) neutrino flux is of vital interest as well as it could reveal 
sources typically associated with very distant and energetic astrophysical objects such 
as AGNs, GRBs and even first signs for new physics beyond the Standard Model (see 
[3J 03 for a review of current models). 

However, the diffuse flux must be observed above the background of atmospheric 
neutrinos. The so called conventional atmospheric neutrino flux from pion and kaon 
decay produced in cosmic ray interactions in the atmosphere falls steeply with energy 

oc Ep, 7 « 3.7 for E v > 10 TeV). At higher energies O(10 5 GeV) the contribution 
of neutrinos from the prompt decay of charmed mesons will begin to dominate the 
flux, thereby hardening its spectral index to 7 ^2.7-3.0. The neutrino contribution 
from charm decay is yet unobserved and predictions vary by more than one order of 
magnitude. A measurement would allow to infer the cross-section for production of 
charmed particles. 

Hence there are a number of interesting potential contributions to the diffuse flux of 
high energy neutrinos. In this paper we address the question under which conditions the 
different contributions can be resolved in future measurements. A Likelihood function 
is constructed incorporating both observable energy and angular information as well 
as systematic uncertainties of the background. The information content of future 
measurements is then analyzed by means of the Fisher matrix technique. We confirm 
the well know fact, that a cubic kilometer sized detector such as IceCube has a large 
discovery potential. However, at the same time we show that it also has a significant 
discrimination potential with respect to the various models. We analyze the information 
content contained in the various detection channels, namely the cascade- and muon- 
channels, and show what can be gained from their combination. 

This paper is organized as follows. In the next section the various fluxes, 
experimental detection channels as well as the simulation of event rates are being 
discussed. In section |3] we review the Fisher matrix method and present its 
implementation with respect to the problem at hand. In section HI we show the results 
of this analysis for a number of different cases. We conclude with a discussion of the 
results and compare to previous work done on similar subjects. 

2. Simulation 

The following fluxes will be considered in various combinations throughout the paper. 
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• An astrophysical flux of neutrinos: dN/dE = a agn x E^ 7agn . The slope 7 agn as 
well as the normalization ct agn are treated as unknown parameters. The assumed 
true value for the slope is 7 agn = 2 while a agn is allowed to vary from 1CT 9 
to 10~ 7 cm"~ 2 s -1 sr -1 GeV -1 . The upper limit is chosen slightly below present 
experimental upper bounds [HI IB1 E] • The Waxman-Bahcall (WB) upper bound 
[TU] (see also [II]), which is of significant astrophysical relevance, corresponds to 
«agn = 4 x 10~ 8 cm _2 s _1 sr _1 GeV~ 1 . Neutrino oscillations will result in approximate 
equipartition of the total neutrino flux among all flavors ^2]- «agn refers to the 
normalization of the flux of one flavor. 

• A flux of neutrinos from charm decay in the atmosphere. This flux is approximated 
by a power-law spectrum with 7 c harm = 2.8. The normalization a c harm ranges from 

to cm -2 s -1 sr -1 GeV -1 , reflecting the range of available models around 
100 TeV neutrino energies jlSl EH E3- Although the shape of the spectrum is 
generally rather well predicted, small variations arise for example from the use of 
different cross-section parameterizations. Hence we again leave both normalization 
and slope as free parameters. A similar flux is assumed for v e and with no 
contribution from v T . 

• The atmospheric neutrino flux. Here we use a selected model [IB] and we do 
not assume any free parameters. It was shown that the uncertainty in the flux 
prediction degrades the ability to detect additional flux contributions [EJ. A 15% 
systematic uncertainty in the normalization is assumed [THj . The normalization is 
allowed to vary as a function of energy with a correlation length of one decade 
in energy. The assumed correlation in energy of the uncertainty reflects the 
fact that relevant ingredients of the calculations, as for example the cosmic ray 
composition or the interaction cross-sections of mesons in the atmosphere, are not 
anticipated to show sudden unaccounted variations. Hence in an experiment, the 
systematic uncertainties between two neighboring energy bins is correlated. The 
implementation of this correlated uncertainty is discussed in sectional 

Note that neutrino oscillations have a very small effect on the atmospheric neutrino 
flux, due to the high energies considered here. 

• The atmospheric muon flux. An analytical parameterization was used [TH] . which 
even at energies of 10 7 GeV agrees within a factor of two with a full CORSIKA 
Monte Carlo simulation [20] • The flux of atmospheric muons at the energies of 
interest (>PeV) is less well predicted, because of the high energies involved. We 
assume an uncertainty of 50 %, again with a correlation length of one decade in 
energy. 

Neutrino events ranging from 10 3 to 10 11 GeV are simulated using the Monte-Carlo 
simulation program ANIS [21] . ANIS allows generating //-events of all flavors, propagates 
them through the Earth and finally simulates //-interactions within a specified volume. 
All relevant Standard Model processes, i.e. charged current (CC) and neutral current 
(NC) i/iV-interactions as well as resonant v e — e~ scattering are implemented. Neutrino 
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regeneration as expected in NC-scattering, v + N — > v + X, and in r production and 
decay chains, v T + N — > r + X, r — ► v T + (z/j) + X, are included at all orders. The 
density profile of the Earth is chosen according to the Preliminary Earth Model [22J. 
Deep inelastic v — iV-cross-sections were calculated in the framework of perturbative 
QCD (pQCD). Parameterization of the structure functions were chosen according to 
CTEQ5 [2H], with logarithmic extrapolations into the small-x region. Tau decay was 
simulated using the TAUOLA program [2*4*] . 

The muon energy after propagating a distance X is approximated by: E^ = 
(E^q - a/b)e~ bX - a/b, with a = 3 MeVg^cm 2 and b = 4 x 1(T 6 g^cm 2 [TJ. Note 
that we are neglecting the fluctuations in muon range leading to slightly over-optimistic 
sensitivities for muons. 

An IceCube-like detector is simulated. It is embedded in ice at a depth of 2 km 
with a bedrock starting in 3 km depth. We define three classes of events, all of which 
have been used by the AMANDA collaboration to search for a diffuse flux of neutrinos: 
cascades, upgoing muons and downgoing muon events. These classes are explained 
below. 

• Cascade events which consist of electro-magnetic or hadronic cascades (also named 
showers). The particle cascade is localized in space and hence can only be observed 
near or within the detector volume. Charged current u e and u T interactions as well 
as NC interactions of any neutrino flavor lead to cascade-like events. In case of NC 
interactions, the visible energy is just a fraction of the incoming neutrino energy. In 
case of a CC tau neutrino interaction the vertex cascades and the cascade obtained 
from tau decay are spatially separated [2H|. It is not yet established under which 
conditions and at which energy this double bang signature can be experimentally 
resolved, hence tau neutrino identification will not be assumed. 

A main advantage of the cascade channel is the reduced background of atmospheric 
muons. An analysis performed with AMANDA has shown that for this class of 
events an almost uniform directional sensitivity can be reached [SI EE!- This leads 
to an improved sensitivity at higher energies, where Earth absorption effects become 
relevant. It is assumed that the background of atmospheric muons, due to their 
different event topology, can be totally eliminated above 1 TeV cascade energies. 

• Upgoing muon events from CC interactions producing an energetic muon 
transversing the detector. These events are selected by a zenith angle cut, 
9 > 90°, to eliminate the background of atmospheric muons. This is the traditional 
observation mode. The vertex of such muon events can be many kilometers away, 
hence the event rate for this class of events is generally higher than in the case 
of cascades. A further advantage is that for these events an angular resolution of 
about 1 degree can be achieved |2]. At high energies, neutrino absorption reduces 
the event rates in this channel significantly. 

• Downgoing muon events from CC interactions do not suffer of absorption in the 
Earth, hence they allow to extend the sensitivity of the instrument to ultra high 




Figure 1. Observable energy spectrum at the detector site. An energy resolution 
of CTiogB = 0.2 was assumed. Left: The cascade energy released by electron neutrino 
interaction within a 1 km 3 sized detector for three different input spectra and averaged 
over all directions. Also shown is the contribution expected from tau neutrinos with 
a flux similar to that of 7 = 2 electron neutrinos. Middle: The upgoing muon event 
rate as a function of the muon energy at the detector. A detection area of 1 km 2 was 
assumed. Right: The downgoing muon event rate as a function of the muon energy 
at the detector. Additionally shown is the event rate due to atmospheric muons. 
The two input power-law spectra have slopes of 7 = 2.0(2.8) and normalizations 
a = 4x 10~ 8 (5 x 10~ 4 )cm- 2 s- 1 sr- 1 GeV- 1 . 



energies beyond 10 9 GeV [BIEHI- The disadvantage of this channel is that there is a 
large background of atmospheric muons. Therefore the effective energy threshold of 
this channel is 10 6 — 10 7 GeV. Here we will assume that the direction of these events 
can be reconstructed with a resolution of better than 5°. This is not yet achieved 
for AMANDA, however the large increase in size and modern readout technology 
should easily allow IceCube to achieve this resolution. 

The resulting event rates are shown as a function of the observable energy in figure^ 
The left plot shows event rates for electron neutrinos for different input spectra. Also 
shown are event rates for tau neutrinos from an astrophysical flux with 7 agn = 2.0. 
The "bump" around 6 x 10 6 GeV is due to resonant u e — e~ scattering. For the highest 
energies the rate of tau neutrinos drops below that of electron neutrinos. That is because 
of the increasing tau decay length. The tau decays more frequently outside ( "behind" ) 
the detector, hence it does not contribute to the visible energy. The middle plot show 
the event rates for upgoing neutrino-induced muons as a function of the muon energy 
in the detector. The right plot show the event rates for downgoing neutrino-induced 
muons as well as for atmospheric muons. 

An energy resolution of a\ Qg E ~ 0.2 for all event classes is assumed. Current 
algorithms used for the operating AMANDA detector are already achieving this 
resolution for cascade- like events [HI EE]- For muon events this is a mild extrapolation 
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of the currently obtained resolution: \ogE^ « 0.3 [213 EDI- However, the results of this 
work are not very sensitive on the exact value of the energy resolution. 

3. The Method 

One of the most popular and powerful methods for parameter estimation is the maximum 
Likelihood method, in which the parameters are determined by maximizing a Likelihood 
function L with respect to its parameters p. Because of the low neutrino event rates, 
the appropriate Likelihood function is based on Poisson statistics: 

M n i e -ni 

Hp) =up(Ni\mm = n^F^ 

where i is the index of the data bin, iVj is the number of observed events in data bin 
i and is the number of events expected given the parameters p. It is convenient to 
define the Log-Likelihood function as C = — log L so that for the Likelihood function 
defined in Equation one obtains: C = —J2i^i^ nn i ~ n i — In AT,!. 

Here we are interested in the achievable accuracy of possible future measurements. 
To evaluate the sensitivity of the upcoming instruments one assumes a fiducial 
model, represented by p, and analyzes the expected averaged constrains on the input 
parameters. However, there is more than one way to do this. One can simulate a future 
experiment in detail and then average over a large number of experimental outcomes by 
means of Monte Carlo technique. Such an approach produces very accurate results but at 
the same time it is computationally demanding and cumbersome in its implementation. 

Here we chose a different method to calculate the sensitivity, namely the Fisher 
matrix method [22j. This method has become a popular tool in modern cosmology 
(see for example j2SJ), most likely due to its transparency and computational efficiency 
Curiously, it is rarely used in the field of astroparticle physics. Hence we briefly review 
this method here (see [2H] for a more in depth discussion of the method). 

The Fisher information matrix is defined as the ensemble averaged Hessian matrix 
of the Log-Likelihood function evaluated at its minimum, 

It represents the information content of the Likelihood function in the close vicinity 
of the true parameters. There are two important inequalities: a Vi < 1/ y/Tu, a statement 
about the lower limit obtainable on the uncertainty a Vi on the parameter pi in the case 
that all other parameters are exactly known. This theorem is known as the Cramer-Rao 



inequality. The second inequality is: o n < y^F -1 )^, which is applicable in the case of 
all parameters being unknown (the parameters j ^ i are marginalized over). Here we 
will assume that a Pi = yj (F~ r )u, which is a reasonable approximation for well behaved 
Likelihood-functions. 

For the Poisson based Likelihood function (flj). one obtains the simple expression: 

F.. = V— — — (3) 

13 * U k d Pi d Pj 



Measuring Diffuse Neutrino Fluxes with IceCube 7 

Note that so far only statistical errors have been assumed. What follows is 
an attempt to incorporate anticipated systematic uncertainties in the background 
prediction. We assume that the conventional atmospheric neutrino flux normalization 
at all energies is known with a precision of cr atmo = 15% [TH], and that this uncertainty 
has a correlation length of one decade in energy. The uncertainty alters the neutrino 
event rate prediction n' { = (1 + r/j)nj, where rji is a nuisance parameter reflecting the 
systematic uncertainty. Its probability function, which needs to be multiplied with the 
right hand side of Equation (0), is assumed to be Gaussian: 

m oc with Cil = 2 ^ (1 + || 1 0g( ^ )|) , (4) 

The dimension of the Fisher matrix increases by the number of nuisance parameters. 
Assigning the indices i, j to the physical parameters of interest (e.g. flux normalization or 
spectral slope) and p, r) to the nuisance parameters, one obtains the following additional 
Fisher matrix elements: 

F m = E C n + S m { ^^- and F lS = £ (5) 
k n k Y n k VP* 



4. The sensitivity 

As described in section El we assume a neutrino flux composed of three components: 

<j>(E„) = <P{E u ) atmo + a am E~^n + a chaim E-^*™ , (6) 

with 7 agn , 7 c harm, ^agn and a c harm being the free parameters of interest. Additionally there 
is a flux of atmospheric muons, which is the relevant background for the downgoing muon 
channel. 

For this analysis the data was binned in energy-bins of the size of the energy 
resolution: Alog-E = 0.2. It was found that for the cascade and upgoing muons 
channel, there is little sensitivity increase when binning the data in energy and zenith 
angle direction. However, in the case of downgoing muons, the sensitivity generally 
increases by about a factor of two, if the reconstructed zenith angle is used as additional 
information (with bins of A cos #=0.1). This is because of the increase in overburden 
towards horizontal directions. The background of atmospheric muons decreases fast 
towards the horizon, while the signal rate increases due to increased target material. 
Hence for the downgoing muon channel explored in this work, the additional directional 
information is being used. 

Figure El shows the significance a/a a = aJ (F -1 )^ with which the normalization 
a can be measured as a function of a assuming three full years of data. The 
various lines correspond to the different detection channels. The dotted line represents 
the combined measurement, where the corresponding Fisher matrix isf ^combined = 

-^cascade ^up-muon ^down-muon' 

| Combining the individual Fisher matrices in a simple sum is only valid if there are no correlations 
between the channels. Certainly there will be some common systematic uncertainties, so the simple 
sum should be considered an optimistic assumption. 



Measuring Diffuse Neutrino Fluxes with IceCube 

100.0 



b 

CD 

u 
c 
o 
u 

c 

00 



cascades 
combined 



muon-up 
muon— down 



10.0 

1 .0 

0.1 
10 









3d- 




1 a 1 




a =5x10" 5 ] 

charm 




10" 
a 



10~1 7 0~ 5 



10" 
a 



10" 



agn 



charm 



Figure 2. Significance (denned as a/a a ) of the measurement of a as a function of 
a. Left: The significance for a agn assuming an additional contribution of neutrinos 
from charm decay at the level of a c harm = 5 x 10 _5 cm _2 s _1 sr _1 GeV _1 . Right: The 
significance of a c h ar m assuming an extraterrestrial contribution with slope 7 agn = 2 
and normalization a agn = 8x 10 _9 cm _2 s _1 sr _1 GeV _1 . 



The left figure shows the significance of the extraterrestrial component with a charm 
contribution of strength a c harm = 5 x 1CT 5 cm _2 s _1 sr _1 GeV _1 while the right figure 
shows the significance of the measurement of the charm component with an assumed 
extraterrestrial contribution of strength a agn = 8 x 10~ 9 cm _2 s _1 sr -1 GeV _1 . 

These curves have a simple interpretation. We assume a measurement is called 
a discovery once the signal is established with 3-cr significance (or in other words 
incompatible with zero at the 99.73 % level). One can read off the corresponding a 
for which an ensemble of experiments will report a discovery. For smaller values of 
a one will typically obtain only upper limits, while for larger values of a a discovery 
becomes more significant. 

As can be seen, a flux at the WB bound would be clearly detected in all three 
channels. The cascade channel, being most sensitive, will allow to make a detection at 
the a = 6 x 10~ 9 level, while the combined sensitivity still improves the results notably. 

A flux of neutrinos from charm decay will be detected in all channels only for the 
highest of all normalizations. Here the cascade channel is significantly more sensitive 
allowing to detect a flux for a c harm a s low as 5 x 10~ 5 cm _2 s _1 sr _1 GeV _1 . However, in 
the presence of a significant flux of astrophysical neutrinos with spectral index 7^2, 
it will become difficult to probe the lowest of all model predictions. 

Note that a discussion on the expected uncertainties on the normalization itself 
becomes more complicated because of the strong covariance between the spectral shape 
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Figure 3. Averaged error ellipse in the 7 agn — 7charm plane for the different channels. 
The columns from left to right show l-,2- and 3-sigma contours. The different 
rows show various combinations of flux normalizations. The diagonal lines represent 

Tagil = '"/charm- 



and the normalization a. The covariance is large, since a is defined for 1 GeV while 
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the relevant energy scale§ E s is 100 TeV to 1 PeV. Redefinition of the normalization 
to a s = aE] would reduce the covariance. By presenting the sensitivity as a/a a , the 
results are independent of the energy scale assumed. 

Once a signal has been established the shape of the underlying spectrum becomes 
of interest. Figure El shows the error ellipses in the 7 agn — 7 c harm plane obtained for 
different combination of small/large extraterrestrial/charm fluxes. The normalization 
values used are a agn = (0.8 — 4) x 10 -8 cm~ 2 s _1 sr _1 GeV _1 and a c h arm = (0.4 — 4) x 
10 -4 cm _2 s _1 sr _1 GeV _1 . The different columns show different significance contours 
ranging from 1- to 3-sigma (from left to right). The diagonal lines represent 7 agn = 
7charm- The interpretation of these ellipses are - within the framework of the Fisher 
matrix method - the following: Given the true values of (a agn , 7 agn ) and (a c harm, 7charm), 
future measurements of these parameters will be distributed according the probability 
contours. 

In case the result of a measurement intersect the diagonal line, separation of the 
two components becomes impossible. As can be seen in the top row of the the figure, if 
both contributions are large, separation is possible at the 3-sigma level. If a agn is large, 
but a c h a rm is small, only 7 agn can be meaningfully constrained. This is not surprising, 
considering that a c h arm has not yet been measured with 3-sigma significance. The other 
two combinations both show at most a 2-sigma separability of the two components. 
Again, of the individual channels it is the cascade channel which provides the tightest 
constrains. 

5. Discussion 

This paper gives an outlook on the sensitivity of future measurements of the diffuse 
neutrino flux. The primary detection channels used by AMANDA, which have 
proven capabilities [HI El U\ , have been investigated and their sensitivity for measuring 
normalization and spectral shape of various potential fluxes was evaluated for a kilometer 
sized detector. 

Thereby it was found that the cascade channel has the highest sensitivity. With 
three years of data, it should be possible to discover a flux of extraterrestrial neutrinos 
at a level as low as 4 x 10~ 9 E~ 2 cm" 2 s _1 sr _1 GeV _1 . It was found that the downgoing 
muon channel is rather sensitive, too, if in addition to muon energy the directional 
information is used. The traditional way to search for high energy neutrinos, namely 
by identifying upgoing muons, is least sensitive. This is because neutrino absorption 
in the Earth suppresses the observable neutrino spectrum at higher energies, therefore 
reducing significantly the lever arm needed to determine the spectral shape. 

A combination of the individual channels leads to the largest sensitivity. 
Nevertheless, for confirmation it would be desirably to detect a significant signal in 
at least two of the three individual channels. 

§ A representative energy scale can be defined as the weighted average of the energy of the bins, with 
the weights given by the bins contribution to the Fisher matrix. 
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A subset of the results presented here has been discussed previously. A detailed 
simulation of the muon channel for IceCube has resulted in precise estimates of the 
sensitivity to a single component diffuse flux |2j. However, spectral shape determination 
and the separability of an astrophysical flux contribution from a contribution from 
atmospheric prompt charm decay were not discussed. It was pointed out before that 
the cascade channel has a significant potential f3J 123 CZ] ■ The ability to determine the 
spectral shape using the cascade channel was studied in Hooper et al. [T7j. Because 
neutrino propagation in the Earth was not included in their simulation the analysis 
was restricted to downgoing cascades only. Further only CC interactions were taken 
into account. Finally the energy bins were chosen to be one order of magnitude, hence 
significantly larger than the energy resolution and bin size assumed here, reducing the 
sensitivity further. This list of differences explains why comparable sensitivities are 
obtained assuming only three years of data here while ten years of data was assumed in 

D3I- 

The work presented here is a first attempt to compare the various detection 
channels in a quantitative manner. Future analyzes might improve upon the sensitivities 
presented here. If tau neutrinos could be identified as such it would provide an important 
additional signature [3T], as there is negligible background due to atmospheric tau 
neutrinos. Hence a single event would be a significant discovery. In practice the effective 
volume/area of an analysis is not confined to the geometrical one. For energetic - hence 
bright - events it could be larger, leading to higher sensitivities. Finally, one should 
restate that the analysis presented here assumes the availability of three years of data, 
while IceCube operation is planned for a time period of 10 years. 
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